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ABSTRACT 



Newly emerging magnetic flux can show a compli- 
cated linked or interwoven topology of the magnetic 
field. The complexity of this linkage or knottedness 
of magnetic flux is related to the free energy stored 
in the magnetic field. Magnetic reconnection pro- 
vides a process to release this energy on the time 
scale of the dynamics. At the same time it approxi- 
mately conserves the total magnetic helicity. There- 
fore the conservation of total magnetic helicity is a 
crucial constraint for the relaxation of complex mag- 
netic fields. However, the total magnetic helicity is 
only the first, most elementary, quantity of an infi- 
nite series of topological invariants of the magnetic 
field. All these invariants are strictly conserved in 
ideal magnetohydrodynamics. As an example a pre- 
liminary set of these invariants is derived. The rel- 
evance of these higher order invariants for the final 
state of relaxation under magnetic reconnection and 
their implications for the release of magnetic energy 
are discussed. 

Key words: magnetic fields; topological invariants; 
magnetic reconnection. 



1. INTRODUCTION 



Magnetic fields in the corona often show a non-trivial 
topology, that is the magnetic field lines are interwo- 
ven or knotted. They form, in mathematical terms 
"knots" and "links". Here "knot" refers to a single 
field line of non-trivial topology, while "link" is used 
if there are at least two field lines which can not be 
separated without cutting of lines. Instead of using 
isolated field lines we can apply these notions to flux 
tubes as well. Note that the existence of complicated 
knots and linkages of magnetic flux is in no way ar- 
tifical. Consider for instance a closed magnetic flux 
tube. In magnetostatic equilibrium the surfaces of 
constant pressure are magnetic surfaces, that is the 
magnetic field is everywhere tangent to the surface. 
Depending on the ratio of poloidal to toroidal compo- 
nent the field lines will close after a certain number of 
windings along the flux tube (n) and around the core 
of the flux tube (m), or they will ergodically fill the 
whole flux surface. In the former case the flux surface 
is called a rational surface, while in the latter case its 
called an irrational surface according to the termi- 



nology in tokamak physics. In generic cases each of 
the infinite set of these rational surfaces has differ- 
ent ratio n/m and therefore for each rational surface 
the field lines are knotted in form of a so called torus 
knot of type (n,m). A simple example of these torus 
knots is shown in Figure \m). 




Figure 1. (a) A thin flux tube forming a torus knot of type 
(n=2,m=5) and (b) a simple linkage of two magnetic flux 
tubes, which leads to a non-vanishing total helicity and 
thus to a lower bound of the magnetic energy. 



Magnetic fields which show knotted or linked mag- 
netic flux contain a certain amount of free energy. 
Here the notion "free energy" denotes the difference 
between the magnetic energy E{B) — J v _B 2 /(87r)<i 3 r 
of initial configuration and the lowest possible mag- 
netic energy, that is the energy of a vacuum magnetic 
field satisfying the same boundary conditions. This 
free energy is stored in the magnetic field and could 
be set free in a relaxation process. The process of 
relaxation can be formally separated into two subse- 
quent processes. First an ideal relaxation, that is a 
relaxation under an ideal Ohm's law, which conserves 
the magnetic topology and leads to a lower magnetic 
energy. Second a non-ideal relaxation, for instance 
by magnetic reconnection, which might finally lead 
to the vacuum field and thus to the lowest energy 
state. Of course in reality most likely the reconnec- 
tion process will set in before the lowest energy state 
accessible under an ideal evolution is reached. More- 
over, it will probably not directly lead to the vac- 
uum field but to a higher energy state, which again 
might be the starting point of a subsequent relaxation 
process and might involve further reconnection pro- 
cesses. But nevertheless the difference between the 
lowest energy state accessible by an ideal evolution 
and the energy of the corresponding vacuum state is 
well determined and is refered to as minimum free 
energy. It is only determined by the topology of the 



initial state and boundary conditions. 

During a non-ideal relaxation the minimum free en- 
ergy can be converted in thermal or kinetic energy 
of the plasma. Thus it might be the source of the 
energy which is required to heat the solar corona and 
the determination of the free energy could give an 
important estimate for the energy available for this 
process. In the next section we will show the relation 
between the free energy and the magnetic topology 
which should motivate the investigation of higher or- 
der topological invariants. A preliminary version of 
such measures of the magnetic topology will be given 
in Section 3, while the third section is devoted to the 
evolution of these invariants under reconnection. 



2. LOWER BOUNDS FOR THE FREE ENERGY 



Roughly speaking the minimum free energy increases 
with the complexity of the magnetic field. This quali- 
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H(B) = / A-Brf 



(1) 



can provide such a lower bound. (For a detailed 
discussion of magnetic helicity see Brown, Canfield, 
|Pcvtsov 1999|) . In the simplest case, that is if the 
magnetic held is enclosed in the volume V (no mag- 
netic flux is penetrating the boundary dV of the vol- 
ume , B • ii\qv = 0), this reads 



E(B) > C\H(B)\, 



(2) 



where C is a constant depending only on the shape 
and size of the volume V. Hence, if the magnetic 
field has a non-vanishing total helicity, resulting for 
instance from a simple linkage of two flux tubes as 
shown in Figure 0b), the magnetic energy can not 
decrease below this bound as long as the topology of 
the field is conserved. A lower energy can only be 
obtained in an evolution which changes the magnetic 
topology, like for instance magnetic reconnection. 
The conservation of topology is usually provided by 
the ideal evolution of the plasma (E + v x B = 0), 
although the condition is more general and there 
are several cases of non-ideal evolution which do not 



chang e the magnetic topology (Hornig & Schindlcr 
1996). 



However, the total helicity is only a very rough mea- 
sure of the topology of a configuration and the lower 
bound given in Eq. |2j might be very low or even zero 
although the topology of the field is non-trivial. An 
example for a configuration which has a vanishing 
total magnetic helicity but is non-trivially linked are 
the so called Borromean rings shown in Figure Ha) 
and in a topological equivalent configuration in Fig- 
ure Bb) . For this configuration the lower bound given 
by Pj) is zero, but a generalization of such a n in- 
equality was given by Freedman fc He (1991b) with 
the help of the so called asymptotic crossing number, 
which shows that there exists a non-zero lower bound 
for the energy in this case as well. 




b) 



Figure 2. (a) The Borromean rings. A configuration of 
vanishing total helicity but non-trivially linked, (b) An 
topological equivalent configuration in form of braided flux 
tubes. 



A drawback of the use of the so called asymptotic 
crossing number is its complicated and abstract def- 
inition. It requires the decomposition of the field in 
closed flux tubes (which is not possible in general) 
and a complicated minimization process applied to 
each combination of linked flux tubes. While this 
number is important to prove the existence of lower 
bounds of the energy, it is impossible to calculate it 
for a generic magnetic field. This is in sharp contrast 
to the total magnetic helicity, which can easily be 
calculated for arbitrary magnetic fields irrespective 
of whether a decomposition into flux tubes exists or 
not. 

This is the motivation to look for invariants, analo- 
gous to magnetic helicity, which give non-vanishing 
values for higher forms of knottedness or linkage 
of magnetic flux in ideal (topology conserving) evo- 
lutions. These higher order topological invariants 
would provide inequalities similar to Eq. |^. An ad- 
ditional reason for the search for such invariants is 
that they would lead to a deeper understanding of 
the global properties of magnetic fields. 



3. HIGHER ORDER TOPOLOGICAL 
INVARIANTS 



A formula for higher order topological invariants 
analogous to the magnetic helicity does not exist yet. 
However, there are constructions of such invariants 
if the magnetic field is confined to a set of isolated 
magnetic fl ux tubes. Suc h a formula was given for 
instanc e in [Berger (19901 for a third ord er invariant 
and in Akhmetiev & Ruzmaikin (1995) for a forth 
order invariant. INote that we are looking for invari- 
ants which depend on the magnetic field only, while 
in the presence of other conserved quantities, such 
as mass density or entropy, it is possible to com bine 
these and in this way construct new invariants (Tur 
k Yanovsky 1993| ). 



In the following a method is provide to construct 
an infinite sequence of higher order invariants which 
can be calculated for arbitrary magnetic fields, but 
have the drawback that they require a specific gauge. 
Therefore, they should be considered as a preliminary 
form, which, however, allows us to gain insight in the 
properties of such invariants. 



The invariants we are looking for are of the form 



H (n) = I h (»)( B ) d 3 
Jv 



(3) 
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for a magnetic field B(x,t) enclosed in the (simply 
connected) volume V. These are topological invari- 
ants if for a frozen-in magnetic field, 



d t B - V x (v x B) = 0, 
the densities h^ n ' satisfy 

d t h {n) + V-(v h {n) ) = 0. 
This results for the integral in 



(4) 
(5) 



gauge for A and G, which for instance does not allow 
to integrate A by 



A = 



V, 



1 



x B(r') d 



3„' 



(17) 



since this corresponds to a different gauge (V • A = 
0), which in general does not satisfy Eq. pi Thus, 

the value of is not uniquely determined, but 

constant for a frozen-in magnetic field. 



(6) 



for a comoving volume V. It is easy to check that 
Eq. H holds for defined by 

:= A • G< n > (7) 

provided the vector fields A(x,t) and G(x,t) satisfy 

d t A + V(v-A) - v x (V x A) = (8) 
9 t G(") - V x (v x GW) + v V'GW = 0. (9) 

Now, if A is the vector potential of B, Eq. [| yields 

d t A - v x V x A = Vx- (10) 



Thus we can meet Eq. || with the gauge A(x, t) 
A(x, t) + V*(ie, t) defined by 



d t ® = -x-v-(A + V#) 



(11) 



Note that this equation can be integrated in time for 
an arbitrary initial gauge "i>(x,0). Similar, defining 

GW by 

V • G (n) := h"- 1 (12) 
Eq. ^| for n — 1 results in 



<9 t G (n) +v V-G (n) = V x F 



(13) 



for some vector field F(x,t). Again we can meet Eq. ^| 

using a gauge G^{x,t) = G^(x,t) + V x J(x,t) 
defining J by 



dtJ := -F + v x (G (n) + V x J), 



(14) 



and the gauge at an initial time is free (J(x, 0)), but 
the equation determines the gauge for all later times. 

We can now rename A by A and G by G. Starting 
with G^ 2 ) = B equations (^) and (O) yield an infi- 
nite recursively defined sequence ofintegrals of the 
magnetic field, all of which satisfy Eq. pf. The first 
(H^) is the total magnetic helicity, which is a sec- 
ond order invariant in that it depends quadratically 
on B 

H {2) = [ A B d 3 r. (15) 
Jv 

The next invariant is of third order 

#(3)= / A-Gd 3 r with V-G (n) :=A-B, (16) 
Jv 

and correspondingly is of order n. The draw- 

back of this construction is the choice of the special 



4. RELAXATION UNDER RECONNECTION 



To release the above defined minimum free energy of 
a magnetic field it is necessary to change its topol- 
ogy. The most important process for astrophysical 
plasmas which allows for a change of the magnetic 
topolo gy is magne tic reconnection. Taylor conjec- 
tured (Taylor 1974) that the total helicity should be 
approximately conserved during a relaxation process 
which involves reconnection. This conjecture proved 
to be true in that the total helicity is decreasi ng on a 
longe r time scale than the magnetic energy (Bergei 
1984). Moreover, there exists a special form of re 



connection events which exactly conserve the helicity. 
This is the case if exactly oppositely directed mag- 
netic fields reconnect such that the magnetic field is 
locally two dimensional and the electric field is per- 
pendicular to the magnetic field. Then E-B = and 
the source term on the right hand side of the balance 
equation for the helicity density, 



dAB 

dt 



6B + E x A) 



2 E B 



(18) 



vanishes. Integrated over a comoving volume with a 
surface everywhere tangential to B and on which the 
evolution is ideal E = — v x B, Eq. |l^ yields 



dt 



A B d A r = -2 



E-B d 3 r. 



(19) 



Thus a reconnection process which satisfies E-B = 
exactly conserves the total helicity. But also for the 
case E • B ^ the total helicity is usually a well 
conserved quantity since reconnection processes in a 
plasma are strongly localized and thus only a small 
fraction of the volume contributes to the integral on 
the right hand side of Eq. || flHornig 1999| ). This 
means that the lowest energy state accessible under 
reconnection is not the vacuum field, but a field which 
has the same value of the total helicity as the initial 
field. This reduces the minimum free energy for this 
kind of relaxation and consequently we have to ask 
whether higher forms of linkage or knottedness will 
lead to additional constraints and thus to a further 
reduction of the minimum free energy. For the case of 
reconnection with E-B = one can easily prove that 
this is not the case. The argument is that for this case 
magnetic reconnection is a two dimensional process 
and can be represe nted by a simple cut and paste of 
magnetic flux (see Freedman & Berger 1993). Using 
this picture of reconnection one can transtorm any 
magnetic field consisting of a set of isolated linked 
or knotted flux tubes into a set of unlinked and un- 
knotted flux tubes which have at most an internal 
twist corresponding to the non-vanishing total helic- 
ity (see Fig. pi). This confirms the original conjecture 



cS 




Figure 3. For the example of Fig. ||| a series of six recon- 
nection processes leads to an unlinked field. 



of Taylor that the total helicity is the only conserved 
quantity. 

Since the higher order linkage or knottedness is not 
conserved for the almost ideal magnetic reconnection 
with E B = 0, we do not expect it to be conserved 
for E-B 7^ 0. However, from the above example 
of (preliminary) topological invariants it is obvious 
that those quantities are integrals of high order of 
the magnetic field. They are global properties of the 
magnetic field and with each integration, that is with 
increasing order of the invariants, more and more in- 
formation of a particular local geometry of the field 
is lost. On the other hand the magnetic reconnection 
process is a local process driven by the local geom- 
etry of the magnetic field. Thus it is questionable 
whether the relaxation process is effective in destroy- 
ing higher forms of knottedness and linkage. With 
other words, to resolve a complicated link or knot re- 
quires a sequence of reconnection processes, and the 
order and locations of these reconnection processes 
is not arbitrary (see Fig || for an example). There- 
fore the relaxation might not lead at the right place 
in the volume to a locally strongly sheared magnetic 
field, which would be necessary to trigger a reconnec- 
tion process and resolve the linkage. Moreover, the 
number of reconnection processes needed to resolve 
a complicated topology increases with complexity. 
Thus these invariants, although they may fluctuate 
under reconnection events, might not decay quickly 
in generic situations. Note that this does not apply to 
situations as in most technical plasmas, where the re- 
laxation process or the non-ideal region, respectively, 
occupies the whole volume under consideration. Only 
if length scales and the magnetic Reynolds numbers 
are as high as in astrophysical plasmas magnetic re- 
connection can be considered as a local process. Also, 
one has to be cautious in applying these considera- 
tions to open field configurations, because in this case 
all forms of linkage and knottedness have to be de- 
fined with respect to a reference field (e.g. the vacuum 
field). 



5. SUMMARY 



A complex magnetic field topology can store large 
amounts of magnetic energy and its minimum mag- 
netic energy under ideal relaxations increases with in- 
creasing complexity Hence such field structures can 
serve as an efficient energy reservoir. However, the 
complexity also imposes restrictions on the release of 
energy from this reservoir. This is because magnetic 
reconnection is the only process capable to release 
this energy on a short time scale. But with respect 
to the first topological invariant, magnetic helicity, 
reconnection turns out to be inefficient to relax cer- 
tain configurations to their lowest energy state (the 
vacuum field), since it approximately conserves this 
quantity. This is not true for higher order topolog- 
ical invariants, as can be shown easily. However, a 
complete relaxation in these cases requires a special 
sequence of reconnection processes, which is unlikely 
to occur in real plasmas. Thus the higher forms of 
linkage do not provide a strict constraint for the final 
state of relaxation by reconnection, as does magnetic 
helicity, but they will significantly delay the process. 
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